Typeset with jpsj2.cls <ver.1 .2> 



Full Paper 



Enhanced Kondo Effect in an Electron System Dynamically Coupled with Local Optical Phonon 

Takashi Hotta 



o 
o 

<N 

Oh 
< 
O 

m 



a 

O 

o 



> 
o 

00 

m 

o 
o 



" 1— "H 



Advanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195 

(Received February 1, 2008) 

We discuss Kondo behavior of a conduction electron system coupled with local optical phonon by 
analyzing the Anderson-Holstein model with the use of a numerical renormalization group (NRG) method. 
There appear three typical regions due to the balance between Coulomb interaction U ee and phonon- 
mediated attraction U p h- For U ee > U p h, we observe the standard Kondo effect concerning spin degree 
of freedom. Since the Coulomb interaction is effectively reduced as U ee — U p h, the Kondo temperature 
Tk is increased when U p h is increased. On the other hand, for U ee < U p h, there occurs the Kondo effect 
concerning charge degree of freedom, since vacant and double occupied states play roles of pseudo- spins. 
Note that in this case, Tk is decreased with the increase of U p h- Namely, Tk should be maximized for 
Uee ~ Uph- Then, we analyze in detail the Kondo behavior at U ee = U p h, which is found to be explained 
by the polaron Anderson model with reduced hybridization of polaron and residual repulsive interaction 
among polarons. By comparing the NRG results of the polaron Anderson model with those of the original 
Anderson-Holstein model, we clarify the Kondo behavior in the competing region of U ee ~ U p h- 

KEYWORDS: Kondo effect, Anderson-Holstein model, Polaron, Numerical renormalization group 
method 



1. Introduction 

Kondo effect and its related phenomena have been cur- 
rently investigated intensively and extensively in the research 
field of condensed matter physics, 1 even after more than forty 
years have passed since the pioneering work of Kondo in 
1964. 2 It has been widely recognized that the Kondo-like phe- 
nomenon generally occurs in a conduction electron system in 
which a localized entity with internal degrees of freedom is 
embedded. Then, a new mechanism of Kondo phenomenon 
with non-magnetic origin has been potentially discussed, al- 
though the original Kondo effect concerning local magnetic 
moment has been perfectly understood. 

Concerning such non-magnetic Kondo effect, Kondo him- 
self has first considered a conduction electron system which 
is coupled with a local double- well potential. 3,4 Two possi- 
bilities for electron position in the double- well potential play 
roles of pseudo-spins and the Kondo-like behavior is consid- 
ered to appear in such a two-level system. In fact, it has been 
shown that the two-level Kondo system exhibits the same be- 
havior as the magnetic Kondo effect. 5 ' 6 Recently, four- and 
six-level Kondo systems have been also analyzed, 7-9 in or- 
der to understand magnetically robust heavy-fermion phe- 
nomenon observed in SmOs4Sbi2. 10, 11 

The multi-level Kondo problem is considered to stem from 
Kondo physics in electron-phonon systems. For instance, the 
present author has discussed how Kondo-like phenomenon 
occurs in a conduction electron system coupled with local 
Jahn-Teller phonon. 12-15 In order to overview the situation, 
it is convenient to envisage the electron potential in an adia- 
batic approximation, although in actuality, the potential is not 
static, but it dynamically changes to follow the electron mo- 
tion. When we simply ignore anharmonicity in the potential of 
Jahn-Teller phonon, there exists continuous degeneracy along 
the circle of the bottom of a Mexican-hat potential, character- 
istic of Jahn-Teller system. When we further take into account 
the effect of cubic anharmonicity, three potential minima ap- 
pear in the bottom of the Mexican-hat potential. Then, we ef- 



fectively obtain the three-level Kondo system. 

It is quite meaningful to pursue a new possibility of Kondo 
effect in the multi-level Kondo system or in the Anderson 
model dynamically coupled with local Jahn-Teller phonon. 
On the other hand, it is also important to confirm the fun- 
damentals of Kondo physics in electron-phonon systems, in 
parallel with the research of complex models with close re- 
lation to actual materials. Thus, we believe that it is useful 
to clarify the Kondo behavior of a simple electron-phonon 
model. In this sense, here we pick up the Anderson model 
coupled with local optical phonon, called the Anderson- 
Holstein model. 16-18 Concerning the origin of heavy-fermion 
phenomenon in SmOs4Sbi2, the periodic Anderson-Holstein 
model has been also analyzed, 19 and a mechanism of the mass 
enhancement due to electron-phonon interaction has been ad- 
dressed. 

We note that the effect of Holstein phonon on the Kondo 
phenomenon is considered to be limited, if we use the adia- 
batic approximation, since the adiabatic potential would sim- 
ply reduce the Coulomb interaction. It is rather interesting 
to examine dynamical phonon effect on the standard Kondo 
behavior concerning spin degree of freedom. From a view- 
point of the relation with actual materials, the Kondo behav- 
ior of the Anderson-Holstein model in the anti-adiabatic re- 
gion may have possible relevance with the enhanced Kondo 
effect in molecular quantum dots. 20-24 Since such a system 
is composed of light atoms, relatively high frequency phonon 
may play an important role for the Kondo physics through the 
competition with Coulomb interaction. 

In this paper, we analyze the Anderson-Holstein model 
by using a numerical renormalization group (NRG) method. 
The results are classified into three categories, depending on 
Coulomb repulsion U ee and phonon-mediated attraction U p h- 
For U ee > E/ph, we easily understand that the standard Kondo 
effect occurs, since the repulsive interaction is still dominant, 
even though it is effectively reduced as U ee — f7 p h. Then, the 
Kondo temperature Tk is increased when /7 p h is increased. 
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Note that in this paper, Tk is defined as a temperature which 
exhibits a peak in the specific heat. For U ee < LT ph , we ob- 
serve the Kondo effect concerning charge degree of freedom, 
since vacant and double occupied states play roles of pseudo- 
spins. In this case, Tk is decreased with the increase of L7 ph . 
Thus, T K is considered to be maximized for U ee « LV ph . 
We focus on the case of U ee = /7 ph and the results are 
found to be understood by the polaron Anderson model with 
reduced hybridization and residual Coulomb interaction be- 
tween polarons. We discuss in detail the Kondo behavior of 
the Anderson-Holstein model in the region of U ee = /7 ph in 
comparison with the polaron Anderson model. 

The organization of this paper is as follows. In Sec. 2, we 
introduce the Anderson-Holstein model and provide a brief 
explanation of the NRG technique used here. We also show 
the canonical transformation of the model to diagonalize the 
phonon part, in order to visualize the competition between 
Coulomb repulsion U ee and phonon-mediated attraction I7 ph . 
In Sec. 3, we show our numerical results for U ee > /7 ph and 
U ee < LT ph . In order to analyze the Kondo temperature, we 
introduce the effective s-d models for both cases. In Sec. 4, 
we discuss in detail the NRG results for U ee = /7 ph . For the 
purpose of intuitive understanding of the results, we propose 
the polaron Anderson model. Finally, in sec. 5, we briefly 
discuss the difference in the Kondo effects between Holstein 
and Jahn-Teller phonons. Throughout this paper, we use such 
units as H=kB=l and the energy unit is set as eV. 

2. Model and Method 

2. 1 Anderson-Holstein model 

Let us introduce the Anderson model coupled with local 
optical phonon. The model Hamiltonian is expressed as 

H = J2 £ kc{ a c ka + ^2(Vcl a d a + h.c.) + H locj (1) 

k<j k<j 

where e k is the dispersion of conduction electron, c ka is an 
annihilation operator of conduction electron with momentum 
k and spin a, d a is an annihilation operator of localized elec- 
tron on an impurity site with spin a, and V is the hybridiza- 
tion between conduction and localized electrons. We choose 
V=0.25 and the energy unit is a half of the conduction band- 
width, D, which is set as 1 eV throughout this paper. 
The local term iTi oc is given by 

#loc = UeeTl^ni + lip + iTeph, (2) 

where U ee denotes Coulomb interaction, n a =d^d a , p is a 
chemical potential, and p=n^ + n^. We adjust p appropriately 
to consider the half-filling case, but the explicit value will be 
shown later. 

The electron-phonon coupling term i^ e ph is given by 

#eph = gQp + P 2 /2 + uj 2 Q 2 /2, (3) 

where g is the electron-phonon coupling constant, Q is nor- 
mal coordinate of breathing mode, and P is the correspond- 
ing canonical momentum. Note that the reduced mass of the 
breathing mode phonon is set as unity. Using the phonon op- 
erator a defined through Q=(a + a))/ \/2cJ, we obtain 

#e P h = wy/a{a + a ] )p + uo(a ] a + 1/2), (4) 

where a is the non-dimensional electron-phonon coupling 
constant, given by a=g 2 /(2u 3 ). The phonon basis is given 



by \t)=(a^Y\0)/V^, where £ is the phonon number and |0) 
is the vacuum state. In actual calculations, the phonon basis is 
truncated at a finite number, which is set as 400 in this paper. 

2.2 Numerical renormalization group method 

In this paper, the Anderson-Holstein model is analyzed by a 
numerical renormalization group (NRG) method, 25 in which 
momentum space is logarithmically discretized to include ef- 
ficiently the conduction electrons near the Fermi energy and 
the conduction electron states are characterized by "shell" la- 
beled by N. The shell of N=0 denotes an impurity site de- 
scribed by the local Hamiltonian. The Hamiltonian is trans- 
formed into the recursion form as 

i^Ar+i = ^AH N + t N y m ]( c N<T C N+i<r + cjv +la CAr a ), (5) 

where A is a parameter for logarithmic discretization, ca/v 
denotes the annihilation operator of conduction electron in the 
TV-shell, and tjq indicates "hopping" of electron between 7V- 
and (N + l)-shells, expressed by 



tN = 



(1 + A- 1 )(1-A- 7V - 1 ) 



2y/(l - A- 2Ar -!)(l - A- 2Ar - 3 ) ' 
The initial term Hq is given by 

H = A~ 1/2 [H loc + J2 V(clda + dUoa)}. 



(6) 



(7) 



The free energy F for local electron in each step is evalu- 
ated by 

F = -T(lnTie- HN/T - lnTre-^ /T ), (8) 

where a temperature T is defined as t=A~^ n ~ 1 ^ 2 in the 
NRG calculation and denotes the Hamiltonian without 
the hybridization term and iTi oc . The entropy is obtained 
by Si mp =—dF/dT and the specific heat Ci mp is evaluated by 
C imp =-Td 2 F/dT 2 . In the NRG calculation, we keep m low- 
energy states for each renormalization step. In this paper, A is 
set as 2.5 and we choose ra=5000, but for some case, it is nec- 
essary to increase m up to 7500 to obtain convergent results. 

In order to clarify the low-temperature properties, we eval- 
uate charge and charge susceptibilities, given by, respectively, 



Xc 



and 



1 ^ e ~E n /T _ e -E m /T 



z ^ 

n,m 



E — E 



\{n\(p-{p))\m)\\ (9) 



Xs 



1 ^ e 

n,m 



-E n /T 



-Em/T 



E — E 



-\(n\a z \m)\ 2 . 



(10) 



where E n is the eigenenergy for the n-th eigenstate \n) of 
H, Z is the partition function given by Z=^ n e~ En ^ T , (p)= 
(l/Z) e~ En l T (n\p\n), and cr z =n^ — n^. We perform the 
calculation in each step by using the renormalized state. 

2.3 Lang-Firsov transformation 

The NRG calculations are performed for the Anderson- 
Holstein model, but in order to grasp roughly the property of 
the model, it is convenient to employ the Lang-Firsov canon- 
ical transformation, 26 defined through the change of an op- 
erator A into A=e R Ae~ R with R= v /ap(a j[ — a). Then, the 
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Fig. 1. (Color online) (a) Txc and Txs vs. temperature for the Anderson Fig. 2. (Color online) (a) T%c and T%s vs. temperature for a=0.3, o;=0.5, 
model with U ee = 1. (b) *Si mp and C- imp vs. temperature for the Anderson and U ee =l. (b) *Si mp and C- imp vs. temperature for the same parameters 



model with U ee = 1. 



as (a). 



Anderson-Holstein model is transformed into 

kcr ka 

+ (LT ee - LT ph )n T n| + (/i - o^)p + u{o) a + 1/2), 

(11) 

where X=e^ a_at ) and /7 ph =2acj. Note that the Coulomb 
repulsion is reduced by the phonon-mediated attractive inter- 
action. We easily understand that for U ee > f7 p h, the local 
state is doubly degenerate with spin degree of freedom, while 
for U ee < U p h, the vacant and double occupied states are 
degenerate at half-filling. From the canonical transformation, 
we obtain that the chemical potential at half-filling is given 
by /i=— /7 ee /2 + 2auo. In the following section, we will dis- 
cuss the numerical results in three regions as U ee > J7 p h, 
U ee < ^ P h, and Uee = U ph . 

3. Kondo Behavior for U ee ^ Z7 p h 

3. 1 Numerical results 

First let us briefly review the results of the Anderson model, 
in order to clarify the effect of Holstein phonon. In Fig. 1(a), 
we show Txc and T\ s for Z7 ee =l. The charge susceptibility 
is rapidly suppressed due to the effect of on-site Coulomb in- 
teraction, while the spin susceptibility is increased. As is well 
known, in the Kondo system, the renormalization flow moves 
toward the strong-coupling regime in which the spin suscep- 



tibility is suppressed, through the local moment regime with 
enhanced spin susceptibility. 

The Kondo behavior is clearly observed in the entropy and 
the specific heat. After the charge susceptibility is suppressed 
around at a temperature in the order of U ee , we can observe 
the local moment region with log 2 between 0.01 < T < 0.1. 
Then, the entropy of log 2 is gradually released and around at 
T ~ 10 -4 , it eventually goes to zero. We can see a clear peak 
in the specific heat around at T ~ 10 -3 due to the release 
of spin entropy log 2. We define the Kondo temperature 
of the Anderson model as a lower peak in the specific heat. 
It is well known that is scaled by e - 1 /( 2 PoJ r o) ? w here the 
exchange interaction Jo is given by J =4V 2 /U ee and po is 
the density of states at the Fermi level. Although a peak in the 
specific heat does not indicate exactly the Kondo temperature, 
the effect of the prefactor is simply ignored here, since in this 
paper, we concentrate on the scaling relation when the pa- 
rameters of the model are changed. Thus, we conventionally 
define Tk as the lower peak in the specific heat throughout 
this paper. 

Now let us show the results of the Anderson-Holstein 
model. First we consider the case in which the Coulomb in- 
teraction is dominant. In Fig. 2(a), we show the results for 
T\ s and T\ c for i7 ee =l, a=03, and co=0.5. We note that 
in this case, Z7 ph =0.3, which is smaller than U ee . We find 
that the charge susceptibility is rapidly suppressed around at 
T ~ 0.03. With decreasing temperature, the spin susceptibil- 
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Fig. 3. (Color online) (a) T% c and Txs vs. temperature for a=l. 5, lj=0.5, 
and U ee =\. (b) Simp and Ci mp vs. temperature for the same parameters 
as (a). 



ity is suppressed around at T ~ 10 -4 . In Fig. 2(b), we show 
the entropy and specific heat for the same parameters as in 
Fig. 2(a). At high temperatures as T > 0.1, we observe the 
entropy larger than log 4 due to the low-lying phonon excita- 
tion states. Passing through the narrow local moment region, 
the entropy log 2 is released and the peak appears in the spe- 
cific heat. Note that the peak position is slightly shifted to the 
higher- temperature side in comparison with Fig. 1(b). 

The overall features of Figs. 2 are quite similar to those of 
Figs. 1 except for the high-temperature region, since the sit- 
uation is effectively understood by the Anderson model with 
reduced Coulomb repulsion U ee — Z7 p h=0.7. As we will dis- 
cuss later, due to the decrease of Coulomb interaction, the 
Kondo temperature is increased in comparison with that of 
the Anderson model with U ee =l. 

Next we consider the situation with U ee < U p h- In 
Fig. 3(a), we show the results for T\ s and T\ c for U ee =l, 
a=1.5, and c^=0.5. In this case, we obtain f7 p h=1.5, leading 
to on-site attractive interaction with U ee — U p h=— 0.5. Thus, 
in sharp contrast to Fig. 2(a), the spin susceptibility is first 
rapidly suppressed around at T ~ 0.1. With decreasing tem- 
perature, through the localized charge state, the charge sus- 
ceptibility vanishes around at T ~ 10 -6 . 

In Fig. 3(b), we depict the entropy and specific heat for 
the same parameters as in Fig. 3(a). At high temperatures, 
we again observe that the entropy is larger than log 4 due to 
the low-lying phonon states. We find the localized charge re- 



gion with log 2 between 10 -4 < T < 0.1, originating from 
the double degeneracy of vacant and double occupied sites. 
With decreasing temperature, the entropy log 2 is eventually 
released and a clear peak appears in the specific heat. 

Due to the effective attractive interaction, we can find the 
Kondo behavior concerning charge degree of freedom for 
U ee < U p h- Since we are now considering the half-filling 
case, as mentioned above, vacant and double occupied states 
are degenerate and these states play roles of pseudo-spins. 
Thus, we can understand the occurrence of the Kondo-like 
behavior in the case of U ee < U p h- The difference in the 
magnitude of the Kondo temperature between the cases of 
U ee > U p h and U ee < U p h will be discussed in the next 
subsection. 

3.2 Effective s-d models 

In the previous subsection, we have shown the NRG results 
for Uee > U p h and U ee < U p h- Here let us discuss analytic 
expressions for Tk. By using the second-order perturbation 
theory in terms of the hybridization, we can derive the ef- 
fective s-d model from the Anderson-Holstein model. There 
appear the virtual second-order processes concerning phonon 
excitations in addition to electronic excitations, which affect 
on the exchange interactions. The similar calculations have 
been done in the derivation of the effective model from the 
Hubbard-Holstein model. 27-29 

After some algebraic calculations, we obtain the effective 
s-d models for U ee > U p h. and U ee < J7 p h, respectively, as 



H s-d = J2 £kC ka C ka 



k,k' 



(12) 



and 



ka 



k,k' 



4i ( 



where S z =(n^ - n x )/2, S+=d\d i , S-=d\d h ^=(|2)(2| 

|0)(0|)/2, 77+=|2)(0|, and r/-=|0)(2| with \2)=d\d\\0). 
The exchange interactions are expressed as 



J ± = W 2 e~ a 



1 



£=0 



t\ \u-u 



ph| 



2£co' 



(14) 



We note that this expression does not hold around at U ee = 
Uph- Note also that the longitudinal and transverse parts of 

(c) (s) 

Hs-d are gi yen by J+ and J_, respectively, while for # s _ d , 
both are given by J + . The difference between J + and J_ 
clearly appears in the asymptotic form for large a as J + ~ 
1/a and J_ ~ e~ 2a /a. Namely, in the strong electron- 
phonon coupling region, J_ decays very rapidly, while J + 
becomes small slowly. The smallness of J_ for large a origi- 
nates from the immobile nature of bi-polaron. Single polaron 
can be relatively mobile in comparison with bi-polaron. 
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symbols denotes numerical results while solid curves indicate analytic re- 
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suits of ! and y . The open circle at ai=l indicates that the analytic 
results cannot be used at this point. 



Now we discuss the Kondo temperatures on the basis of the 

(s) 

effective s-d models. For the case ofC/ee > E7p h , H^ A is the 
isotropic s-d model and we can easily obtain 

T K )= M-2^)- (15) 

(c) 

On the other hand, i? s _ d is the s-d model with anisotropic ex- 
change interaction. For this case, Shiba has obtained the ex- 
plicit expression for the binding energy E. 30 When we define 
the Kondo temperature Tk as Tk = — E, we obtain 



.(c) 



Dexp 



-1 



ztanh 



2poJJl 



J\-Jl 



j2 V J+ 

(16) 

In Fig. 4, we depict Tk/T-^ vs. a. Numerical results are 
shown by solid symbols. Note that Tk as well as in the 
numerical results is defined as a temperature which shows the 
peak in the specific heat. Analytic results for the s-d models 
are depicted by solid curves, which indicate eqs. (15) and (16) 
We find that the numerical results agree well 



divided by t£ 0) 



(s) (c) 

with the analytic curves for H^_ d and H^_ d . Note that the 

numerical results seem to scatter due to the effect of discrete 

temperature A - ^ -1 ^ 2 , but all the results are considered to 

be in the error-bars. 

As mentioned above, for U ee < U p h, the effective s-d 

model becomes highly anisotropic for large a. Since the trans- 
fer 

verse part is exponentially small, the interaction part of if s _ d 
becomes Ising-like and the Kondo temperature is rapidly sup- 
pressed. Thus, the plot of Tk vs. a is asymmetric at the center 
of a=l (U ee =U p h)- Note also that two analytic curves seem 
to converge to a value similar to the numerical result at a=l, 
but as mentioned above, the expression of Tk does not hold 
around at U ee =U p h- The Kondo behavior in the repulsion- 
attraction competing region with U ee « U p h will be sepa- 
rately discussed in the next section. 

4. Kondo Behavior in the Competing Region 

Now we focus on the case of U ee = J7 p h. In Fig. 5(a), we 
show the results for T\ s and T\ c for a=l, U ee =l, and co=0.5. 
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Fig. 5. (Color online) (a) T\c and T\s vs. temperature for a=l and o;=0.5 
with U ee =2a(ju. (b) *Si mp and Ci mp vs. temperature for the same parame- 
ters as (a). 



As naively expected from the cancellation of the on-site in- 
teraction, the situation is considered to be understood by the 
non-interacting model. In fact, as shown in Fig. 5(a), it is very 
difficult to distinguish the charge and spin susceptibilities, in- 
dicating that the interaction effect is virtually ignored. 

In Fig. 5(b), the entropy and specific heat are shown for 
the same parameters as in Fig. 5(a). At high temperatures 
as T > 0.1, in common with the cases of U ee ^ J7 p h, we 
find that the entropy is larger than log 4 due to the low-lying 
phonon states. We note that the entropy gradually goes to zero 
without showing the region of log 2. We can see the peak in 
the specific heat around at T ~ 0.01, which has been assigned 
as the Kondo temperature in Fig. 4. However, this Kondo- 
like behavior is not due to Coulomb interaction, but it origi- 
nates from the hybridization in the non-interacting Anderson 
model. This point will be discussed later again. 

Next we increase the value of a by keeping the relation of 
U ee = U p h- In Fig. 6(a), we show the results for T% s and 
T\c for a=9, Z7 ee =2acc;, and co=0.5. From the viewpoint of 
the comparison with actual materials, the value of a=9 seems 
to be unrealistically large, but we consider such a situation in 
order to complete the discussion from the theoretical view- 
point. As naively expected from the cancellation of on-site 
interaction, the situation can be understood from the non- 
interacting model, but in the strong-coupling case, the results 
clearly show the interaction effect. For T > 10 -3 , T% c and 
Txs takes a constant value of 0.5. This region can be inter- 
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Fig. 6. (Color online) (a) T% c and T\s vs. temperature for a=9 and uj=0.5 
with U ee =2au>. (b) S[ mp and Ci mp vs. temperature for the same parame- 
ters as (a). Note that in this case, we increase m up to 7500, since conver- 
gence becomes worse due to the strong-coupling nature. 



preted as the free orbital regime, described by the conduction 
electron systems. 25 Around at T ~ 10 -3 , spin and charge re- 
sponse begin to be separated from each other. The charge sus- 
ceptibility is rapidly suppressed, while the spin susceptibil- 
ity grows up, suggesting the local moment regime. After that, 
the spin susceptibility is gradually suppressed and it eventu- 
ally goes to zero around at T ~ 10 -9 , entering the strong- 
coupling regime. 

The above behavior can be clearly found in the entropy and 
the specific heat, as shown in Fig. 6(b). The high-temperature 
behavior is almost the same as that in Fig. 5(b). For 0.1 < 
T < 10 -4 , we clearly observe the region of log 4 due to the 
four-fold degeneracy of vacant, spin-up, spin-down, and dou- 
ble occupied states, corresponding to the free-orbital regime. 
Then, the entropy of log 2 relevant to charge degree of free- 
dom is released and a peak is formed in the specific heat. Fi- 
nally, the residual entropy of log 2 for spin degree of freedom 
is released and we find another peak in the specific heat. 

From the above numerical results, we deduce that the case 
of U ee = U p h is described by the Anderson model with effec- 
tive on-site repulsive interaction. From H, eq. (11), we pro- 
pose the polaronic Anderson model, given by 



kcr 



(17) 
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Fig. 7. (Color online) (a) (n^n^) vs. temperature for several values of a 
with uj=0.5 and U ee =2auj. Inset shows (n^n^) at the lowest temperature 
in the present NRG calculations, (b) A e ff and U e ^ vs. a for u;=0.5 and 

U e e=2aUJ. 



where V e f£ is the hybridization for polarons, U e s is the resid- 
ual Coulomb interaction among polarons, and /x e ff is the 
chemical potential for polarons, given by /i e ff=— U e ^/2. 

By taking the average over the zero-phonon state, we ob- 
tain V; ff as V e ^=V(0\X\0)=Ve~ a/2 . However, we cannot hit 
upon an idea to derive the analytic form of /7 eff . Then, we 
resort to a numerical method to determine f7 e ff so as to repro- 
duce the average value of double occupancy {n^n\) of H. 

In Fig. 7(a), we show (n^n\) for several values of a with 
U ee =2auo. In the inset, we depict (n^n\) at the lowest tem- 
perature which we can reach in the present NRG calculation. 
Irrespective of the vales of a, in the high-temperature region 
as T > 10 -3 , {n^ni) takes the value near the non-interacting 
one, 0.25. When we further decrease the temperature, for 
a < 6, {n^rii) is not largely suppressed and it keeps the value 
about 0.2 even at low temperatures. However, for a > 6, it 
is rapidly suppressed and takes the value much smaller than 
0.25. In fact, as shown in the inset, the low-temperature value 
of {n^ni) is rapidly suppressed around at a ~ 6. 

Readers may feel it strange that the double occupancy is so 
suppressed, in spite of the fact that the on-site interaction ex- 
actly vanishes at U ee = U p h- As mentioned in the explanation 
of the effective s-d model, bi-polaron has immobile nature in 
comparison with single polaron for large a. Due to the differ- 
ence in the mobility between polaron and bi-polaron, single 
polaron state has the energy gain of the exchange interaction, 
while bi-polaron cannot. Thus, the single polaron state is fa- 
vored in the strong-coupling region, leading to the effective 
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occurrence of the residual repulsion between polarons. 

By changing the values of f7 e ff for each values of a, we 
repeat the NRG calculations for H e f?, until we can repro- 
duce the low-temperature value of {n^n\) of H. The results 
are shown by solid symbols in Fig. 7(b). In order to under- 
stand intuitively the extent of correlation effect, we also de- 
pict the line of 7rA e ff, where A e ff is the width of the virtual 
bound state, expressed as A e ^=7rpoV^=7rpoV 2 e~ a . Note 
that the expansion parameter of the Anderson model is given 
by Z7 e ff/(7rA e ff). 31 ~ 34 For a < 6.5, we observe that J7 e ff is 
smaller than 7rA e ff, suggesting that the peak in the specific 
heat should be determined by A e ff . For large a, on the other 
hand, /7 eff is still small, but A e ff is exponentially reduced, as 
observed in Fig. 7(b). Then, in this region, we expect that the 
correlation effect becomes significant. 

In Fig. 8, we summarize the numerical results for the peak 
temperature in the specific heat with the fitting curve de- 
duced from the polaron Anderson model. 35 For a < 6.5, 
as described above, the NRG results are simply scaled by 
A e ff . For a > 6.5, as typically found in Figs. 6, we obtain 
clear two peaks in the specific heat. The higher peak con- 
cerning charge degree of freedom should be characterized by 
[/efr, while the lower one indicates the characteristic tempera- 
ture of the standard Kondo effect, given by e -1 / ( 2p0 Jeff ) with 
J e ff =4V r J r /LT e flF • By adjusting appropriately numerical prefac- 
tors at a=7, we can actually fit the NRG results by /7 e ff and 
e -i/(2poJ r eff) w e believe that the success of this fitting sug- 
gests the effectiveness of the polaron Anderson model, which 
describes the low-energy states of the Anderson-Holstein 
model in the competing region of U ee « U p h- 

5. Discussion and Summary 

We have clarified the Kondo behavior of the Anderson- 
Holstein model by using the NRG method. It has been clearly 
shown that the results are categorized into three classes la- 
beled by U ee > U ph , U ee < U ph , and U ee ~ ^ P h- For 
U ee > U p h, the standard magnetic Kondo phenomenon oc- 
curs with the reduced Coulomb interaction, suggesting the 
increase of the Kondo temperature in comparison with that 
of the Anderson model without the electron-phonon interac- 
tion. For U ee < E/ph, me charge Kondo effect occurs and the 
Kondo temperature is decreased with the increase of electron- 
phonon coupling constant, since the relevant exchange inter- 
action is decreased when we increase a. Note that the effec- 



tive s-d model for this case becomes highly anisotropic, since 
the transverse part is related to the bi-polaron motion with im- 
mobile nature in comparison with single polaron. 

Around at U ee ~ U p h, the Kondo temperature is maxi- 
mized, but for small a, the characteristic energy is simply 
given by the width of the virtual bound state of the non- 
interacting Anderson model. In this sense, we should not 
call the peak in the specific heat as the Kondo temperature 
for small a, but the Kondo-like behavior occurs at relatively 
high temperature, when the strong electron-phonon interac- 
tion competes with Coulomb interaction. We may consider 
possible relevance of this scenario to the enhanced Kondo 
temperature observed in molecular quantum dots. For larger 
a, the effect of residual polaron repulsion becomes signif- 
icant, since the polaron hybridization is exponentially sup- 
pressed. Then, we have found the Kondo singlet formation, 
through the free-orbital and local moment regimes. 

Finally, we provide a comment on the different effect of 
the phonon mode on the Kondo temperature. In this paper, we 
have concentrated on the Holstein phonon and stressed the en- 
hancement of the Kondo temperature. However, this enhance- 
ment depends on the feature of relevant phonon. For instance, 
for the case of Jahn-Teller phonon, the Kondo temperature 
is monotonically decreased with the increase of the electron- 
phonon coupling constant. 1213 For Jahn-Teller phonon, the 
local electron-phonon state has double degeneracy originat- 
ing from clockwise and anti-clockwise rotational phonon 
mode. Such geometrical degree of freedom is characterized 
by J=±l/2, where J is total angular moment, composed of 
electron orbital and phonon angular moments. The conduc- 
tion electrons must screen phonon angular moment in addi- 
tion to electron orbital moment to form the singlet ground 
state with J=0. Thus, the Kondo temperature is decreased 
when we increase a. 

In summary, we have discussed the Kondo effect of the 
Anderson-Holstein model. We have observed that the Kondo 
behavior can be explained by the isotropic s-d model for 
U ee > U p h, the polaron Anderson model for U ee ~ U p h, 
and the anisotropic s-d model for U ee < U p h- The Kondo be- 
havior has been found to be enhanced when U competes with 
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